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We present the first evidence of the decay B° — > p°tv°, using f40 fb _1 of data collected at the 
T(4S) resonance with the Belle detector at the KEKB asymmetric e + e~ collider. We detect a 
signal with a significance of 3.5 standard deviations, and measure the branching fraction to be 
B (B° -> p°7v°) = (5.1 ± 1.6(stat) ± 0.9(syst)) x 10" 6 . 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 



Recent measurements of the CP violating parameter 
sin2fa P, have confirmed the Kobayashi-Maskawa 
mechanism |3| as the origin of CP violation within the 
Standard Model (SM). It is now essential to test the SM 
via measurements of other CP violating parameters. Of 
particular importance are the other two angles of the Uni- 
tarity Triangle, fa and fa. Measurements of fa typically 
rely on time-dependent studies of decays of B mesons 
to light mesons, such as B° — > 7r + 7r~ and p^it^ . Al- 
though these analyses are complicated by the presence 
of penguin amplitudes, isospin analyses can be used to 
extract fa 0]. Recent evidence for direct CP violation in 
B° — > ir + ir~ [j| indicates a sizeable penguin contribution; 
furthermore measurements of the B° —> 7r°7r° branching 
fraction at a level higher than most theoretical expecta- 
tions 0] suggest that much larger data samples will be 
needed for a model-independent extraction of fa from 
the 7T7T system using an isospin analysis. 

Measurements of fa from the pir system rely on knowl- 
edge of the branching fraction of B° — ► p°ir°. The isospin 
analysis depends on this information, along with the CP 
asymmetry, since all the other pir final states have been 
observed 0, @ ■ An alternative technique to extract fa 
uses an amplitude analysis of B° — > 7r + 7r~7r° 9]. Since 
B° — > p°tt° results in this final state, it is essential to un- 
derstand its contribution, as well as possible effects from 
scalar resonances, e.g. an , and nonresonant sources |lCj . 

Recent theoretical predictions for the branching frac- 
tion of B° — » p°tt° are typically around or below 
10~ 6 yjj, while the most restrictive experimental up- 
per limit, recently set by the BaBar Collaboration, is 
B (B° -> p a ir a ) < 2.9 x 1CT 6 at the 90% confidence 
level. In this letter, we present the first evidence for 
B° p°ir°. 

The analysis is based on a 140 fb _1 data sample 



containing 152 x 10 6 B meson pairs collected with the 
Belle detector at the KEKB asymmetric-energy e + e~ 
collider 0|. KEKB operates at the T(4S r ) resonance 
(y^s = 10.58 GcV) with a peak luminosity that exceeds 
1.2 x 10 34 cm~ 2 s _1 . The production rates of B + B~ and 
B°B° pairs are assumed to be equal. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a three-layer silicon vertex de- 
tector (SVD), a 50-layer central drift chamber (CDC), 
an array of aerogel threshold Cerenkov counters (ACC), 
a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a su- 
perconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect Kl mesons and to iden- 
tify muons (KLM). The detector is described in detail 
elsewhere jl '1| . 

Charged tracks are required to originate from the inter- 
action point and have transverse momenta greater than 
100 MeV/c. To identify tracks as charged pions, we 
combine specific ionisation measurements from the CDC, 
pulse height information from the ACC and timing infor- 
mation from the TOF into pion/kaon likelihood variables 
C-k/k- We then require C v / {L^ + Lk) > 0.6, which pro- 
vides a pion selection efficiency of 93% while keeping the 
kaon misidentification probability below 10%. Addition- 
ally, we reject tracks that are consistent with an electron 
hypothesis. 

Neutral pion candidates are reconstructed from photon 
pairs with invariant masses in the range 0.115 GeV/c 2 < 
to 77 < 0.154 GeV/c 2 , corresponding to a window of ±3a 
about the nominal ir° mass, where a is the experimental 
resolution for the most energetic ir° candidates. Photon 
candidates are selected with a minimum energy require- 
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ment of 50 MeV in the barrel region of the ECL, defined 
as 32° < 7 < 129° and 100 MeV in the endcap regions, 



defined as 17° < 7 < 32° and 129° < 7 < 150°, where 
7 denotes the polar angle of the photon with respect to 
the beam line. The ir° candidates are required to have 
transverse momenta greater than 100 MeV/c in the lab- 
oratory frame. In addition, we make a loose requirement 
on the goodness of fit of a 7r° mass-constrained fit of 77 
(xlo)- 

Possible contributions to the 7r + 7r~7r° final state from 



backgrounds are explicitly vetoed for 
+ 



Pi 



charmed (6 — » 

the decays B° — > D~n + , D°tt° and J/tfjir , based on 
the two-particle invariant masses. From Monte Carlo 
(MC) simulation, we find a small combinatorial back- 
ground from b — > c remains. 

B candidates are selected using two kinematic vari- 
ables: the beam-constrained mass M^ c = E^ cam 
and the energy difference AE = Eb — E\ xam . Here, Eb 
and pb are the reconstructed energy and momentum of 
the B candidate in the centre of mass (CM) frame, and 
-Ebcam is the beam energy in the CM frame. We con- 
sider candidate events in the region —0.20 GeV < AE < 
0.40 GeV and 5.23 GeV/c 2 < M hc < 5.30 GeV/c 2 . 
With these boundaries 30% of events have more than 
one candidate, and that with the smallest Xvtx + X^-o is 
selected, where Xvtx i s the goodness of fit of a vertex- 
constrained fit of 7r + 7r~. We define signal regions in 
AE and M hc as -0.135 GeV < AE < 0.082 GeV and 
5.269 GeV/c 2 < A/ bc < 5.290 GeV/c 2 respectively. 

To select p°ir° from the three-body tt + tt~~tt° candi- 
dates, we require the invariant mass to be in the 
range 0.50 GeV/c 2 < m n+n - < 1.10 GeV/c 2 and the p° 
helicity angle to satisfy |cos0£ cl | > 0.5, where 0£ cl is de- 
fined as the angle between the negative pion direction and 
the opposite of the B direction in the p rest frame [l5j . 
Contributions from B° — > p ± 7r =F are explicitly vetoed by 
rejecting candidates with tt^tt invariant masses that fall 
into the window 0.50 GeV/c 2 < m^i^o < 1.10 GeV/c 2 . 
This requirement also vetoes the region of the Dalitz plot 
where the interference between p resonances is strongest. 

The dominant background to B° — > 7r + 7r~7r° comes 
from continuum events, e + e~ — > qq (q = u, d, s, c). Since 
these tend to be jet-like, whilst BB events tend to be 
spherical, we use event shape variables to discriminate 
between the two. We combine five modified Fox- Wolfram 
moments [l^ | into a Fisher discriminant and tune the co- 
efficients to maximise the separation between signal and 
continuum events. We define 9b as the angle of the recon- 
structed B candidate with respect to the beam direction 
in the CM frame. Signal events have a distribution pro- 
portional to sin 2 8b , whilst continuum events are flatly 
distributed in cos0_b- We combine the output of the 
Fisher discriminant with cos 9b into signal/background 
likelihood variables, C s /b, and define the likelihood ratio 
1Z = C s / (C s + Cb). In order to maximise the separation 
between signal and background, we make use of the addi- 



tional discriminatory power provided by the flavour tag- 
ging algorithm developed for time-dependent analyses at 
Belle |l7| . We utilise the parameter r, which takes values 
between and 1 and can be used as a measure of the con- 
fidence that the remaining particles in the event (other 
than 7T + 7r~7r°) originate from a flavour specific B meson 
decay. Events with a high value of r are considered well- 
tagged and hence are unlikely to have originated from 
continuum processes. Moreover, we find that there is no 
strong correlation with any of the topological variables 
used above to separate signal from continuum. 

We use a continuum suppression requirement on r and 
1Z that maximises the value of N s /\/N s + Nb, where N s 
and Nb are the numbers of signal and background events 
contained in the intersection of the AE and M^ c sig- 
nal areas. To obtain N s we use a large statistics sam- 
ple of p°7r° MC, and assume a branching fraction for 
B° — > p°7T of 1 x 10~ 6 . We estimate Nb from a contin- 
uum dominated sideband in data, defined as the union 
of the two regions -0.20 GeV < AE < 0.40 GeV and 
5.23 GeV/c 2 < M hc < 5.26 GeV/c 2 , and 0.20 GeV < 
AE < 0.40 GeV and 5.26 GeV/c 2 < M hc < 5.30 GeV/c 2 . 
We use an iterative procedure to find the optimal con- 
tiguous area in r-TZ space. This method is found to be 
robust against statistical fluctuations in the samples used 
to obtain N s and Nb- The result of the procedure is that 
we select events that satisfy either 1Z > 0.92 and r > 0.70 
or 1Z > 0.35 and r > 0.95, as shown in Fig.Q] In addition 
to optimising N s / \/N s + Nb, this requirement is found to 
improve N s /Nb by a factor of 76. 
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FIG. 1: Distributions of signal (MC) and continuum (side- 
band data) events in r-TZ space. The marked region indicates 
the selection requirement obtained from the optimisation pro- 
cedure described in the text. 

Following all the selection criteria described above, 
the signal efficiency measured in MC is found to be 
(1.91 ± 0.01)%, and we find 73 candidates remain in the 
data, as shown in Fig. Ufa). We obtain the signal yield 
using an unbinncd maximum-likelihood fit to the AE- 
Mbc distribution of the selected candidate events. The 
fitting function contains components for the signal, con- 
tinuum background, b — > c background and the charm- 



less B decays B + 



P + P°, 



B^ 



p + n° and £ H 



7T 7T 
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The possible contribution from other charmless B decays 
is found to be small (0.7 events) using a large MC sam- 
ple (lif. and is taken into account in the systematic error. 
The probability density functions (PDFs) for the signal 
and charmless B backgrounds are taken from smoothed 
two dimensional histograms obtained from large MC 
samples. For B + — > p + p° our MC assumes 100% longi- 
tudinal polarisation 19] . For the signal PDF, small cor- 
rections to MC peak positions (< 0.5 MeV) and widths 
(< 16%) are applied. These factors are derived from 
control samples (B° -> D*~ p+ with D*~ -»■ £j 7r _ , 
D° -> K+-K-, p+ -> tt+tt and B+ -> D°p+ with 
D° — > K + t:^ 1 p + — > tt + tt°), in which we require that 
the 7T° momentum be greater than 1.8 GeV/c in order to 
mimic the high momentum tt° in our signal. 

The two-dimensional PDF for the continuum back- 
ground is described as the product of a first-order polyno- 
mial in AE with an ARGUS function [2(J in M hc . Con- 
tributions from b — > c are also parametrised as a prod- 
uct of two one-dimensional PDFs. Using MC we find 
the AE distribution of this background in the fitting re- 
gion is modeled accurately by an exponential function; 
the Mbc distribution is modeled by the ARGUS func- 
tion. All of the shape parameters describing the con- 
tinuum and b — > c backgrounds are free parameters in 
the fit. The normalisations of B+ -> p+n" (2.0 ± 0.5 
events) and B + — ► tt+tt (2.3 ± 0.5 events) are fixed in 
the fit according to previous measurements 0, 0] , while 
the normalisations of all other components are allowed 
to float. 

The fit result is shown in Fig.^b) and (c). The signal 
yield is found to be 15.1 ±4.8 with 3.6<r significance. The 
significance is defined as y/— 21n(£o/'Cmax), where £ max 
(Co) denotes the likelihood with the signal yield at its 
nominal value (fixed to zero). The backgrounds from b — > 
c and from B + — > p + p° form a peak in the low AE region. 
The fit results for these background sources are consistent 
with the MC expectation, which for B + — > p + p° is based 
on our branching fraction measurement |19|. 

In order to check that the signal candidates originate 
from B° —>■ p°ir° decays, we change the criteria on m n + v - 
and cos 0£ cl in turn, and repeat fits to the AE-Mb c distri- 
bution. The yields obtained in each m v + v - and cos6^ cl 
bin are shown in Fig. |2Id) and (e) . 

We use the cos 0^ e ^ distribution to limit contributions 
from (T7T°, /o(980)tt° and 7T+7r~Tr° (nonresonant), which 
are expected to have similar shapes in this variable. 
We perform a \ 2 fit including components for pseu- 
doscalar — > pseudoscalar vector (PV ~ cos 2 0£ c i)> an< ^ 
pseudoscalar — * pseudoscalar scalar (PS ~ flat) decays, 
for which the shapes are obtained from our p°7r° signal 
MC, and a sample of cttt MC ji^], respectively. We find 
the PS level is consistent with zero, and assign a system- 
atic error due to the possible contribution in our signal 
region of ^5%. The m n + w - distribution supports the 
conclusion that our signal is due to B° — > p°7r°. 



To extract the branching fraction, we measure the re- 
construction efficiency from MC and correct for discrep- 
ancies between data and MC in the pion identification 
and continuum suppression requirements. The correction 
factor due to pion identification (0.89) is obtained in bins 
of track momentum and polar angle from an inclusive D* 
control sample (D*~ -> D°tt-, D° -> K + w~). The re- 
sulting systematic error is ±3.3%. For the continuum 
suppression requirement on r and 1Z, we use a control 
sample B° -> D~ p + with D~ -> K + tt~tt~, p + -> tt+tt , 
which has the necessary feature of being a neutral B de- 
cay to ensure the r behaviour is the same as that of our 
signal. A correction factor of 1.15 is obtained; the statis- 
tical error of this control sample accounts for the largest 
contribution to the systematic error, ±11%. 

We further calculate systematic errors from the follow- 
ing sources: PDF shapes 5% (by varying parameters 
by ±ler); tt° reconstruction efficiency ±3.5% (by com- 
paring the yields of rj — > Tr Tr°Tr and r\ — ► 77 between 
data and MC); track finding efficiency ±2.4% (from a 
study of partially reconstructed D* decays). We use our 
calibration control samples to study possible effects on 
the efficiency due to the AE > —0.2 GeV requirement 
and assign a 2% systematic error. The total systematic 
error due to possible charmless B decays not otherwise 
included is ±5.3%. We repeat the fit after changing the 
normalisation of the fixed B decay components according 
to the error in their branching fractions, and obtain sys- 
tematic errors from the change in the result of ±1%. In 
the case that the normalisations of B backgrounds fixed 
in the fit are simultaneously increased by ler, the statisti- 
cal significance decreases from 3.6a to 3.5er; we interpret 
the latter value as the significance of our result. Finally, 
we estimate the systematic uncertainty due to possible 
interference with B° — > p+Tr^ by varying the m 7r ± 7r a veto 
requirement. We find the largest change in the result (by 
9.3%) when this requirement is removed, and assign this 
as the error. The total systematic error is ±17%, and we 
measure the branching fraction of B° — * p°Tr to be 

B (B° -> p°tt°) = (5.1 ± 1.6(stat) ± 0.9(syst)) x 10" 6 . 

In order to test the robustness of this result, a number 
of cross-checks are performed. We vary the selection on r 
and 1Z. We try numerous combinations of requirements, 
with efficiencies that vary between 1.60% and 2.70%. In 
all cases consistent central values of the branching frac- 
tion are obtained. We also repeat the analysis using a 
looser requirement on the lower bound of AE and obtain 
consistent results. Finally, we select p+Tr^ candidates 
from the Tr + Tr _ Tr° phase-space using the same continuum 
suppression requirement, and measure a branching frac- 
tion for B° — > p+tt^ that is consistent with previous 
measurements [7j. 

In summary, we observe the first evidence, with 3.5cr 
significance, for B° — > p°ir° with a branching fraction 
higher than most predictions |ll| . and a central value 
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FIG. 2: (a) Scatter plot of AE versus Mt, c for the selected candidate events; the box indicates the intersection of AE and Mbc 
signal regions, (b), (c) Distribution of AE(M\> C ) in the signal region of M^ C {AE). Projection of the fit result is shown as the 
solid curve; the dashed line represents the signal component; the dot-dashed curve represents the contribution from continuum 
events, and the dotted curve represents the composite of continuum and B-related backgrounds, (d), (e) Distributions of fit 
yields in m n + n - and cos6^ cl variables for p°7r° candidate events. Points with error bars represent data fit results, and the 
histograms show signal MC expectation; the selection requirements described in the text are shown as dashed lines. 



above the upper limit recently set by the BaBar collab- 
oration 8]. Our result may indicate that some contribu- 
tion to the amplitude is larger than expected, which may 
complicate the extraction of (f>2 from the pir system. 
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